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ABSTRACT 

We present the results of a systematic search for large-scale spectroscopic 
variability in apparently single Wolf-Rayet stars brighter than f ~ 12.5. In 
this first paper we characterize the various forms of variability detected and 
distinguish several separate groups. For each star in our sample, we obtained 
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4-5 high-resolution spectra with a signal-to-noise ratio ~ 100. Our ultimate goal 
is to identify new candidates presenting variability that potentially comes from 
Co-rotating Interaction Regions (CIR). 

Out of a sample of 25 stars, 10 were found to display large-scale changes of 
which 4 are of CIR-type (WR 1, WR 115 WR 120 and WR 134). The star WR 134 
was already known to show such changes from previous studies. Three WN8 
stars present a different type of large-scale variability and we believe deserve a 
group of their own. Also, all three WC9d stars in our sample present large-scale 
variability, but it remains to be checked if these are binaries, as many dust- 
making WR stars are double. Finally, of the remaining stars, 10 were found to 
show small-amplitude spectral changes which we attribute to normal line-profile 
variability due to inhomogeneities in the wind, and 5 were found to show no 
spectral variability, as far as can be concluded from the data in hand. 

Follow-up studies are required to identify potential periods for our candidates 
showing CIR-type changes and eventually estimate a rotation rate for these WR 
stars. 

Subject headings: stars: Wolf-Rayet — stars: variables — stars: winds, outflows 
— stars: mass-loss 



1. Rotation and Spectroscopic Variability of Single Wolf-Rayet Stars 



Stellar ro tation has now been fu l ly incorporated into evolutiona ry models of mas- 
sive stars (eg. iMeynet fc Maeden l2000l : iHeger. Langer fc Woosleyl l2000l ) . prompted by sev- 
eral observational facts that were unexplained by p revious models without rotation (eg. 
Meynet et al.lll994l : IWooslev. Langer fc Wea.verlll993l ) and that pointed to a solution related 
to mixing. For example, iHerrero et al.l (Il992l ) had found that many O stars, and in particular 
fast rotators, show an enhanced He abundance; they called this the helium discrepancy. The 
models including rotation have been very successful in reproducing the observed abundances, 
predicting higher He and N abundances for higher masses and higher rotation velocities. 
Evolutionary models with rotation have also provided a solution for many other previously 
unexplai ned observations such as the well-known blue-to- red supergiant ratio (B/R) prob- 
lem (e.g. iLanger fc Maedeii Il995l ). Indeed, when neglecting rotation, evolutionary models 
were unable to reproduce the rapid increase in B /R with increasing metallicity. As the new 
models predict many more red supergiants at low metallicity coupled with a shorter blue 
supergiant phase, they provided a natural solution. 



It is now clear that rotation is an important factor in our understanding of massive stars 
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and their evolution. It is therefore interesting to compare theoretical predictions with the 
available observational data. For O stars, the latest homogenous, lar ge-scale surveys of pro - 
jected rotational velocities have been carried out by lPennyl (Il996l ) and lHowarth et al.l (119971 ). 
who respectively studied 177 and 373 stars, the latter also including early-type B stars. These 
studies were both based on high-resolution lUE spectra and used a cross-correlation tech- 
nique with a narrow-line template spectrum to measure the rotational speeds. Both studies 
note the absence of very narrow lines for supergiants while some smaller velocities are found 
for main-sequence stars. Instead of interpreting this as a possible increase in rotational ve- 
locities as the star evolves, these authors fav our an additional broadening mechanism such as 
macro-turbulence. In fact, iHowarth I (120041 ) strengthens this assertion by pointing out that 
the shape of the absorption lines expected from pure rotation is a poor match to the observed 
profiles and that a line dominated by transverse turbulence is a much better match. Ignor- 
ing this possible effect (most likely small) of macro-turbulence on the projected rotational 
velocities, their distribution peaks at ~ 80 kms~^ and has a very extended high- velocity tail 
reaching to 400 kms~^. These values are not incompatible with model predictions; on the 
main sequence, the models show that even if the star starts out with a relatively flat internal 
rotation profile, the rotation becomes increasingly differential (but never more than a factor 
of two in velocity) because the convective core is contracting and therefore spinning up, while 
the outer layers are expanding and therefore slowing down. In spite of this, there is a global 
decrease in the rotation rate as a function of time as angular momentum is lost through 
mass-loss. So, depending on their a ge and t heir i nitial rotation rate, the observed velocities 
are quite acceptable. Interestingly, iPennyl (119961 ) notes that, in fact, the higher rotational 
velocities are generally as sociated with lower-mass stars. This is an important result. Indeed 
Meynet fc Maeden (120001 ) predict that the rotation rates decrease more rapidly for higher 
initial rotation rates and for higher mass and therefore higher mass-loss rate stars, in view 
of the increased loss of angular momentum. This agrees very nicely with the observations. 

For WR stars, very slow average surface rotation velocities in the range 20-70 kms~^ are 
predicted for the He-burning phase, depending on the initial rotation velo city on the main 
sequence and on the evolutionary state of the star (IMeynet fc Maederll2003l ). The reason for 
these low values and the narrowness of this range is simple; the amount of mass that must be 
removed from massive O stars to turn them into WR stars is so big that most of the angular 
momentum is lost in the process. However, the differential rotation that steadily increased 
on the main sequence suddenly becomes very pronounced, as the helium core contracts at 
the end of core hydrogen burning. The star ends up with a core spinning faster than the 
initial ZAMS rotation velocity. 



A fast-spinning core, at least in some WR stars, would support the presently favoured 
model to explain long-soft Gamma- Ray Bursts (GRBs). Indeed, the so-called collapsar model 
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f lWooslevI I1993I : iMacFadven fc WooslevI Il999l : iMacFadven. Wooslev fc Hegerl |200l|) involves 
the collapse of a rapidly rotating massive star in its pre-supernova phase. In this model, the 
central core becomes a black hole surrounded by a disk. The accretion of m atter through the 
disk c reate s powerful jets in which the matte r travels at very high speeds. iHeger fc Langer 
( I2OOOI ) and iHirschi. Meynet fc Maederl (120041 ). have followed the evolution of a rotating He 
core to determine if the very rapid core velocities required by the collapsar model could be 
reached in spite of the many modes of transport and redistribution of the angular momentum 
taking place (mass-loss, shear, turbulence, convection, etc). They found that indeed this is 
possible for initial rotation speeds on the main sequence that are typical of that of O stars. 

Ho wever, addin g magn etic fields to the models has brought new elements to the problem. 
Indeed, iHeger et al.l (|2005r) h ave shown that including magnetic torques caused by a Spruit- 
Tayler dynamo (jSpruit 2002) dramaticaly reduces the angular momentum of the core leading 
to periods in the resulting pulsar of 15 ms at birth, which is compatible with the average 
value for common pulsars but much too small to make a long-soft gamma-ray burst (GRB). 

Nevertheless, for a rotating star, efficient chemical mixing can occur and if the mixing 



( Maeder 


1987; 


Laneer 


1992) 



which are required for the development of the Spruit-Tayler dynamo (ISpruitll2002l ). This lack 



of a core-envelope structure prevents the magnetic co re-envelope coupling des cribed above 



and the spin-down of the star's core can be avoided. IWoosley &: Hegerl (120061 ) have shown 
that for a rapidly-rotating star on the main sequence and if the mass-loss in the WR phase 
is reduced by a factor of ten, massive-star evolution can indeed lead to a core with sufficient 
angular momentum to produce a long-soft GRB. The necessar y reduced mass-loss in the WR 



phase favors low metallicity environments for GRBs. Finally lYoon et al.l (120061 ) calculated 
evolutionary grids for massive stars including rotation and the effects of a Spruit-Tayler 
dynamo for different metalicities. Still adopting WR mass-loss rates reduced by a factor of 
ten, they predict that the production of long GRBs should be limited to metalicities smaller 
than Z ~ 0.004. This translates into a peak at a redshift of z=4. 

No reliable direct observation of WR rotation rates has ever been obtained. Attempts 
have been made to estimate them based on the w idth of (what the authors had assurn ed to 
be) photospheric absorption lines in WR spectra (lMasseylll980l : iMassey fc Contilll98ll ). but 
these have been seriously questioned, as wind expansion and turbulence have not been taken 
into account. 

It may be virtually impossible to obtain the rotationnal velocity of WR stars directly 
through photospheric line-profiles but it could still be possible to deduce it from periodic vari- 
ability in photometry or of spectral lines formed in the winds. Indeed, the winds of massive 
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stars have been demonstrated to be highly variable. In particular, one type of structure in the 
wind that accounts for some of these variations, Co-rotating Interaction Regions (CIRs), are 
thought to be closely linked to the rotation of the star. Indeed, these structures are caused 
by perturbations at the base of the wind that propagate through it and are carried by ro- 
tation. This generates a spiral-like struct ure in the density distribution, analogous to those 
commonly observed in the solar corona (IHundhausenl Il972l ) . which ver y likely leads to the 
ubiqu itous Discrete Absorption Components in the spectra of O stars (ICranmer fc Owocki 
19961 ). These are rarely observed in WR P C ygni profiles as their a bsorption components are 
usually saturated (one exception is WR24; iPrinja &: Smith! Il992l ). However, it is believed 
that they also lead t o a very characteristi c , larg e-scale, periodic variability pattern in WR- 
wind emission lines. iDessart fc Chesneaul (120021 ) show the unambiguous S-shape pattern in 
time-series of model spectra of WR emission lines produced as a consequence of these varia- 
tions. Observational support has been obtained for t his type of variability in two clear cases 



through repeated spectrosco pic observations. WR6 (ISt-Louis et al. 
P=3.76 days) and WR134 jMcCandliss et al.lll994j : iMorel et al 



1995 



1999 



Morel et al. 1997 



P=2.25 days) show 

periodic variations without any indication of a companion. The correspondence between the 
data and models is striking when one considers that optical depth effects, which are certainly 
occurring in the observed data, have not been included in the models. This brings strong 
support to the idea that C IRs do in fact t race the ro tation rate of the star at the base 
of the wind. According to iHamann et al.l (jl988l . Il995l ). these stars both have i?* ~ SRq. 
This radius corresponds by definition to a Rosseland optical depth of 20. Assuming that 
CIRs form at this radius, the equatorial rotation speed (at i?*) is ~ 40 km s~^ for WR6 and 
~ 70 km s~^ for WR134, not too different from values predicted by evolutionary models. 
The periods and variability patterns were found to be the same for lines formed at various 
distances in the wind, indicating that the CIRs are governed by angular momentum con- 
servation which results in a coherent structure in the wind. The period therefore provides 
a direct measurement of the rotation rate of the underlying star. To calculate the rotation 
velocity we need to know at what radius the perturbation leading to the CIR originates. If 
this radius is known, we can then calculate the rotation velocity at that radius. Much more 
observations and analysis of variability associated with CIRs are required to reach this goal. 
Therefore, the above velocities are only given as rough estimates for illustrative purposes. 

In an attempt to identify more WR stars displaying CIR- type variability, we have carried 
out a survey of all Galactic WR stars brighter than v ~ 12.5. Our comple te sample includes 



39 southern stars, which will be discussed in IChene fc St-LouisI (l2009al ). and 25 northern 
stars, which we discuss here. The goal of this paper is to present the different types of 
spectroscopic variability found in our magnitude-limited sample of stars in a qualitative but 
also quantitative manner. Our aim is to use these characteristics to establish some guidelines 
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to help identify the spectroscopic variations associated with CIRs. Once identified, these 
stars can then be observed more intensively and, if present, periods can be determined. In 
our second paper, we will build on these basis and use the complete sample to determine 
the frequency of the various types of line-profile variability. Follow-up papers on individual 
stars that have been more intensively monitored are also planned. 

Our northern sample includes all Galactic WR stars visible from the Mont Megantic 
Observatory that are not confirmed as SBl or SB2 binaries. Our targets are listed in Ta- 
ble 1 where the stars' name, spectral type, RA, DEC and v are provided. When avail- 
able, we also provide terrninal velocities measured either from lUE absorption troughs 
( jPrinja. Barlow fc Howarthlll99oh. from near- infrar ed line profiles (iHowarth fc Schmutall992l ) 



or from optical line profiles (IHamann et al.ll2006l ). The last column indicates each star's 



variability status (see below). Note that although some stars in our sample have published 
spectral types that indicate a possible companion (not indicated in the table), none have con- 
firmed orbits. In the case of WR4 and WR 128, only photometric variations have been found 
without radial velocity confirmation of binarity. In the case of WR131, WR143, WR156 
and WR 158, the claim of binarity is based solely on diluted e mission lines (d.e.l.) a nd /or 



the presence of absorption lines in the spectrum (+a). However, iHamann fc Grafenerl (120071 ) 
argue that binary status established from the d.e.l or +a criteria is often misleading as WR 
stars can produce a wide variety of spectra depending on their physical parameters. Those 
authors have demonstrated that in many cases these types of spectra can be reproduced 
entirely with a single-star model. In particular, WR 108, WR 156 and WR 158 were among 
the stars that they reproduced with a single-star model. Our final sample consists of 16 WN 
stars and 9 WC stars. We have included WR 134 for comparison even though it is already 
known to show variability thought to be associated with CIRs. However, we have excluded 
WR 157 because it is too close on the sky to another star that is only ~1 magnitude fainter. 
The separation between the two stars is ~2.2", i.e. a value comparable with the average 
seeing at the mont Megantic Observatory, where our data were obtained. It was therefore 
impossible to isolate the WR star in our 1.5" spectrograph slit. 



2. Observations 

Our observations were obtained at the 1.6m telescope of the mont Megantic Observatory 
in Quebec, Canada during 6 runs in 2001 and 2002. Details are given in Table 2 in which 
the run number, the dates of observations and the Julian Dates (JD) are listed. We used a 
830.8 1/mm grating giving a resolution of AA=1.6 A (3 pix.) combined with a 2048x4096 
EEV CCD. 



Table 1. Our Sample of Single WR Stars in the Northern Hemisphere 



Name Spectral Type^ RA (2000) DEC (2000) v Doo(kms"^) Variability 



WRl 


WN4 


00 43 28.4 


+64 45 35.4 


10.51 


2100'^ 


1900* 




LSV 


WR2 


WN2 


01 


05 


23.03 


+60 25 18.9 


11.33 


3200^ 


1800* 




NV 


WR3 


WN3 


01 


38 


55.63 


+58 09 22.7 


10.70 




2700* 




ssv 


WR4 


WC5 


02 


41 


11.68 


+56 43 49.7 


10.53 








NV 


WR5 


WC6 


02 


52 


11.66 


+56 56 07.1 


11.02 








NV 


WR106 


WC9d 


18 


04 


43.66 


-21 09 30.7 


12.33 








LSV 


WR108 


WN9h 


18 


05 


25.74 


-23 00 20.3 


10.16 


1220^ 


1170* 




SSV 


WRllO 


WN5-6 


18 


07 


56.96 


-19 23 56.8 


10.30 


3190^ 


2300* 




ssv 


WRl 11 


WC5 


18 


08 


28.47 


-21 15 11.2 


8.23 


2415^ 


2300=^ 




NV 


WR115 


WN6 


18 


25 


30.01 


-14 38 40.9 


12.32 


1480^ 


1280* 




LSV 


WR119 


WC9d 


18 


39 


17.91 


-10 05 31.1 


12.41 








LSV 


WR120 


WN7 


18 


41 


00.88 


-04 26 14.3 


12.30 


12253 


1225* 




LSV 


WR121 


WC9d 


18 


44 


13.15 


-03 47 57.8 


12.41 


1300^ 






LSV 


WR123 


WN8 


19 


03 


59.02 


-04 19 01.9 


11.26 


970* 






LSV 


WR124 


WN8h 


19 


11 


30.88 


+16 51 38.2 


11.58 


710" 






LSV 


WR128 


WN4(h) 


19 


48 


32.20 


+18 12 03.7 


10.54 


2270^ 


2050* 




SSV 


WR131 


WN7h 


20 


00 


19.12 


+33 15 51.1 


12.36 


1400-* 






ssv 


WR134 


WN6 


20 


10 


14.20 


+36 10 35.1 


8.23 


1905^ 


1960=^ 


1700* 


LSV 


WR135 


WC8 


20 


11 


53.53 


+36 11 50.6 


8.36 


1405^ 


1500=* 




ssv 


WR136 


WN6(h) 


20 


12 


06.55 


+38 21 17.8 


7.65 


I6O52 


1760=* 


1600* 


ssv 


WR143 


WC4 


20 


28 


22.68 


+38 37 18.9 


11.95 








ssv 


WR152 


WN3(h) 


22 


16 


24.05 


+55 37 37.2 


11.67 


1800^ 


2000* 




ssv 


WR154 


WC6 


22 


27 


17.82 


+56 15 11.8 


11.54 


2700^ 






NV 


WR156 


WN8h 


23 


00 


10.13 


+60 55 38.4 


11.09 


660* 






LSV 


WR158 


WN7h 


23 


43 


30.60 


+61 55 48.1 


11.46 


900* 






ssv 



^Spectral types are from Ivan der Hucht ] ll200lh 
^Prinia. Barlow &: Howartbl l ll990h 
feowarth fc Schmutj lll992l ') 
feamann et al.l ll2006ll 
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Table 2: Observing Runs 



Observing 


Dates 


Julian Date 


run No. 








1 


25/06/2001-05/07/2001 


2452086- 


-2452096 


2 


26/07/2001-02/08/2001 


2452117- 


-2452124 


3 


01/10/2001-07/10/2001 


2452184- 


-2452190 


4 


19/06/2002-24/06/2001 


2452445- 


-2452450 


5 


17/07/2002-26/07/2002 


2452473- 


-2452482 


6 


16/10/2002-24/10/2002 


2452564- 


-2452572 



Our goal is to identify new large-scale variability within our sample of stars. The extra 
features we are looking for on top of the main emission lines are quite large and therefore 
easy to find. Based on our previous experience with WR6 and WR134, we have made the 
assumption that a small number of good quality spectra, well sampled in time are sufficient 
to identify new candidates showing variability possibly associated with CIRs. Of course, to 
quantify the changes and determine the period, a more thorough follow-up will subsequently 
be required. Therefore, we set out to obtain 4-5 spectra per star in our sample. 

The spectra were reduced using the IRAF software package in the standard way. First a 
bias was subtracted, then each image was divided by a flat fleld. After extracting the spectra, 
the wavelength calibration was carried out using a He-Ar calibration lamp. Although our 
original spectral range was 1500 A wide, we were forced to reject the reddest part of the 
spectrum due to a lack of usable lines in the calibration spectra. The flnal wavelength interval 
is AA=4500-5200 A. The signal-to-noise ratio of the spectra ranges from 100 to 120. 

As our spectra are not photometrically calibrated, the flnal step of our data reduction 
procedure was to rectify the spectra. To do so, we have fltted a low order Legendre poly- 
nomial to wavelength regions with no strong spectral lines and divided our spectra by the 
fltted curves. 



3. Variability Search 

The flnal spectra are shown in the top panels of Figures [1] to [13] for the two strongest 
lines present in our wavelength interval. For WN stars these are HellA4686 and HellA4860 
while for WC stars they are ClllA4650 and ClvA5016. For each star we plot the differences 
between individual spectra and the global mean (which is shown in the second panel from 
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o) 



the top). The differences have been shifted vertically for clarity. The Heliocentric Julian 
Date (HJD) of each observation is indicated on the left-hand side of the plot and the scale 
factor of the ordinate is indicated in the top right-hand corner of each panel. The graphs 
are organised by spectral type starting with WN stars from early to late types (Figures [T] 
to [8]) followed by WC stars from early to late types (Figures [9] to [T3!). Note that for WR4, 
WRlll, WR5 and WR154, we do not show the plot of the ClllA4650 line because in those 
cases the line was on the extreme edge of the wavelength interval and therefore there was 
not enough continuum on one side of the line to allow for an accurate enough spectrum 
rectification. This could lead to a false detection of line variability. We prefer to adopt a 
conservative approach in this respect. 
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Fig. 1. — a) top: Montage of the HellA4686 (left) and HelA4860 (right) residuals (individual 
spectra — mean) for WR 2 (WN2). For both cases, the scale factor of the ordinate is indicated 
in the top right-hand corner of the plot. HJD - HJD(2000) is indicated in the y-axis. second 
from top: Mean spectrum, second from bottom: cr-spectrum. bottom: S (99%) spectrum, b) 
The same as a) for WR 3 (WN3ha). 





Fig. 3.— Same as FiglDfor a) WR 128 (WN4(h)) and b) WR 110 (WN5-6). 





Fig. 5.— Same as FiglDfor a) WR 136 (WN6) and b) WR 120 (WN7). 




Fig. 6.— Same as FiglDfor a) WR 131 (WN7h) and WR 158 (WN7h). 




Fig. 7.— Same as FiglUfor a) WR 123 (WN8) and b) WR 124 (WNSh). 
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Fig. 8.— Same as FiglDfor a) WR 156 (WN8h) and b) WR 108 (WN9h). 



o 930.76 




b) 



1500 1500 
velocity (km s"') 



500 1500 
velocity (km s"') 



o -025.76 



WR4 
CIVX501 6 



-1500 1500 
velocity (km s"^) 



Fig. 9.— Same as FigHlfor a) WR 143 (WC4) and b) WR 4 (WC5) for the ClllA4650 (left) 
and ClvA5016 (right) hne-profiles. Note that for WR 4, the montage for the ClllA4650 hne 
is not shown as the data are unusable. 
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Fig. 10.— Same as Figj9]for a) WR 111 (WC5) and b) WR 5 (WC6). Note that the montage 
for the ClllA4650 hne is not shown as the data are unusable. 
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Fig. 11.— Same as Figl for a) WR 154 (WC6) and b) WR 135 (WC8). Note that for 
WR 154, the montage for the ClllA4650 hne is not shown as the data are unusable. 
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Fig. 13.— Same as Figlfor WR 121 (WC9d). 
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In order to search for significant spec t rosco pic variability among our targets, we have 



adopted the formahsm of iFuUerton et al.l (Il996l ). First, we have calculated the temporal 



variance spectrum (TVS), which is a statistically rigorous means of identifying the variable 
parts of our spectra. For spectra, the TVS at wavelength j is defined by: 

where the weights of the individual spectra, Wi, are given by Wi = {-^Y with 

'l^ 



„ 2 I ^ \ „ -2 



a standardized dispersion. Here, Sj is the weighted mean spectrum, aic is the noise in the 
continuum and aij are wavelength correction factors given by = (^)^, (Jij being the 
noise associated with the flux Sij. 

Then, we have established the level of significant variability assuming that our data are 
governed by a reduced chi-squared distribution with N— 1 degrees of freedom. Using the 
standardized dipersion, ctq, for scaling, we therefore adopt a threshold of : croXAr-i(99%). 
Finally, we used Sj(99%) = ^/ i ''^"^'^qV') quantify the level of variability at each wave- 



length: a spectrum that reaches a level of 1, is variable at a 1 a level, if it reaches 3 it is 
variable at a 3 o" level, etc. We have chosen a confidence level of 99%, which means that if 
the value of Sj(99%) reaches 3, we are 99% confident that the spectrum is variable at a 3 cr 
level. We note that as the standardized dispersion, ctq is calculated using the noise level in 
the continuum regions and as the spectra are rectified, there is no real information on the 
variability of the continuum regions. In the bottom panels of Figures [T] to [T^ we have plotted 
the S(99%) spectrum for each star for easy comparison with the spectra plotted above. 

Inspection of Figures [1] to [13] reveals that the majority of the stars in our sample show 
significant variability in at least one of the two main spectral lines in their spectrum. One 
notable exception is the WN2 star WR2. For WR4, WR5, WRlll and WR154, we did 
not detect variability in the weak Civ line, but we only have limited information at our 
disposal since we do not have usable data in the strong Cm line. It is clear that the level 
of variability is quite different from one star to another. The variations range from a few as 
to close to 60(TS in WR115. As the quality of our spectra is relatively uniform, we can, in 
principle, compare the a from star to star, but its value only provides information on the 
significance of the changes with respect to the noise level of our data. 

As the line strengths are very different from star to star, comparing the significance level 
of the variations does not provide useful physical information. A quantity that is potentially 
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more interesting is the fraction of the hne flux that is variable. To characterize this, it is 
best to use a quantity that is representative of the variance in pixel j compared to the noise 
at that wavelength (o".y) rather than use a w eight th a t is lin ked to the noise in the associated 



continuum, as the TVS does. To do this, IChena (120071 ) has defined a modified TVS as 
follows: 

N / _ \ 2 

(2) 

where aoj is the reciprocal of the rms of the signal-to-noise ratio at pixel j of the spectral 
line series. 

In order to quantify what fraction of the line flux (without the contiuum) is variable, we 
divide this value by the line flux, [Sj — !)• Of course this quantity has no meaning outside 
emission lines and in particular in continuum regions. To avoid any problems, we therefore 
limited our caculations to a width corresponding to a relative intensity level of 1.3 of each 
line. Also, for very weak lines (< 1.3 in relative intensity), our ratio becomes meaningless 
because the line flux, the denominator, is too close to zero and the variability level, the 
numerator, is dominated by the noise. In such cases, we have left the panel blank. Finally, 
absorption components of P Cygni profiles are also excluded since in that case the value of 
{Sj — 1) will be negative and cannot be compared directly with what is observed in emission 
lines. 

1/2 ~ 

The values of Oj = {TVSmod)/ /{Sj — 1) are plotted in the middle panels of Figures [T] 
to [T3l Several interesting conclusions can be drawn from these plots. First, for all types of 
stars, both lines show roughly the same level of variability. However, one has to note that 
in the case of WN stars, we compared lines of the same ion and therefore it is likely that 
in a given star, these two lines are formed in similar regions of the wind. The line widths 
are indeed, in most cases, very similar. Another interesting feature of the variability is that 
the calculated values of S are much stronger in P Cygni absorption components than in the 
corresponding emission components, i.e. the variations are more significant. This is perhaps 
not surprising since that part of the profile comes from a much smaller volume of the wind 
than the emission and therefore suffers less from cancelling effects from large-scale changes 
in density or ionisation arising in different parts of the wind. This effect is clearly present in 
the P Cygni profile of WR123 (HellA4686) or WR119 (ClllA4650), for example. 

After a careful examination of Figures [T] to [T31 we divide our sample in three main cate- 
gories: stars showing no profile variability (NV), stars showing small-scale profile variability 
(SSV — cr < 5%) and stars showing large-scale variability (LSV — a > 5%). Note that 
the stars designated as NV are not necessarily stars that do not show spectral variability, 
but rather that no significant variability can be detected at a level of confidence of 99% on 
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spectra having signal-to-noise ratio ~100. Our conclusions are listed in the last column of 
Table 1. The detected subpeaks have widths that range from ~140 kms~^ for the smallest 
isolated structures to ~350 kms~^ for the largest. Unfortunately, our hmited dataset does 
not allow us to search for correlations between these widths and, for example, the level of 
variability or the width of the global line. In the following sections we will discuss in turn 
the variability we found in WN and WC stars. 



3.1. WN Stars 



Significant variability has been found in all but one WN star in our samp le, i.e. WR2, 
the W N star of the earliest type. However, this seems to be an unusual WR star. iHamann et al 
( 120061 ) note that this star has unique line-profile shapes. These authors were unable to re- 
produce its spectrum using their model atmospheres without convolving the normal wind 
lines with a rotation profile with a velocity of 1900 km s~^, which is very close to the star's 
critical velocity. This is, of course, a simplistic way of including the effects of rotation on 
the line profile and is likely to lead to improbable conclusions but it reflects the fact that the 
line-formation process for this star is not totally understood and that the star is unusual in 
that respect. 

For WR3, WR152, WR128, WR136, WR131, WR158 and WR108, variability has 
indeed been detected but at a relatively low level of a ~ 2 — 3% of the line flux. The nature 
of the changes is also very different from the CIR-type variability we first set-out to identify. 
The structures are small and are most likely stochastic changes associated with inhomo- 
geneities in the win d of WR stars, such as d iscussed from an observational point of view by 
Moffat et al.l (119881 ) and iLepine et al.l (119961 ). for example. More recently, there have been 
some i nsightful theoretica l investi ga tions i nto srn all-scale inhomogeneities in the winds of hot 
stars. IPessart fc Owockil J2002al fl l2003l . boOSal lbh present radiation hydro dynamic simula- 



tions of hot-star winds line-profile variations that arise as a consequence of the intrinsically 
unstable nature of the line-driving mechanism of the wind. Those calculations show, for 
example, that including Rayleigh- Taylor or thin-shell instablities together with some lateral 
averaging of the diffuse radiation is important to understand the behaviour of the width of 
the subpeaks as a function of their position on the line profile, i.e. the lateral versus radial 
size of the small scale inhomogeneities in the wind. We believe that the characteristics of the 
variations seen in the stars that we classify as SSV, are compatible with the type of changes 
caused by such small-scale inhomegeneities in the wind. 



For WRl, WRllO, WR115, WR134, WR120, WR123, WR124 and WR156, the 
variations have a much larger amplitude. Therefore, it is quite likely that the physical 
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phenomenon from which they originate is different from what is observed for SSV stars. We 
discuss these stars in more detail below. 

WR134 is well- known to show CIR-typ e variability for wh ich a period of 2.3 days has 



been identified (e.g. iMcCandliss et al.l 1 1994 iMorel et al.l Il999l ). We use it as a test that 



indeed with 4-5 spectra we can identify easily the type of variability we are looking for. One 
can clearly see on top of both the HellA4686 and HellA4860 lines, broad emission excesses 
at different positions on the line. The changes are highly significant and reach a = 6-8 % of 
the line flux. 

Three other WN stars show a similar type of variability: WRl, WR115 and WR120. 
These stars clearly show CIR-type variations; large-scale structures are seen superposed 
on the broad emission lines reaching a ~7-8% of the line flux in WRl and WR120 and 
a ~10— 12% of the line flux in WR115. We therefore classify these stars as LSV and further 
put them in the CIR-type category. 

The situation is less clear for WRllO. Clear structures can be distinguished on both 
spectral lines. The variations reach a level of a ~5% of the line flux, which is significantly 
higher than the other stars we have classified as SSV and that we associate with small-scale 
wind inhomogeneities. However, the detected structures are not quite as strong as for those 
showing variability thought to be associated with CIRs. Therefore, we classify this star as a 
SSV and use it to establish our approximate upper limit for the level of small-scale /blob-like 
variability. 

Finally, for WR123, WR124 and WR156, all WN8-type stars, a different type of vari- 
ability is found. The changes reach a relatively high level (from 5 to 10% of the line flux) 
but appear to be different in nature from those thought to be associated with CIRs. The a 
profile has a very typical shape with two maxima, one on the red and the other on the blue 
side of the emission line and a minimum at line center. The only exception is the HellA4686 
profile of WR123 but in that case there is a very strong contrast between the variability in 
the absorption component and that of the emission. In fact this effect of a more variable 
absorption component seems to be amplified in WR108, the only WN9 star in our sample. 
Indeed, despite the fact that we have classified this star as SSV, the changes are different in 
nature from that observed in the other stars we have put in that category. In this case, the 
variability is actually confined to the P Cygni absorption component of the HellA4860 line, 
which is clearly present. In view of the specific nature of the changes, we create a separate 
class of LS V for WN8-type st a rs. Ac tually, WN8 stars are well-known to display an unusual 



behaviour. iMarchenko et al.l (119981 ) presented a comprehensive study of the photometric 
and spectroscopic variability patterns of 10 Galactic WN8 stars and found a link between 
the level of variability in the continuum and in spectral lines which they suggest may be 
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lin ked to pulsationa l insta bilities. Other unusual behaviours of WN8 stars as summarized 
by lMarchenko et al.l (119981 ) include a low binary frequency, a high distance from the Galactic 
plane and/or runaway speeds and the fact that they are rarely found in stellar clusters and 
associations. Therefore, it is perhaps not surprising to find that they fall in a class of their 
own compared with other types of WR stars when it comes to spectral variability. 



3.2. WC Stars 



Because of problems with our spectrum rectification procedure of the ClllA4650 lines 
in our early-type WC stars, the information at our disposal for these stars is more limited 
than for WN stars. Indeed, this line is located very close to the blue extremity of our 
wavelength interval. Therefore, the continuum regions that are available bluewards of the 
line are much more restricted. This reduces the accuracy that can be achieved during the 
spectrum rectification procedure. Of course, the higher the contrast of the line with respect 
to the continuum, the more it is affected by this effect. The WC4 star WR143 has a 
somewhat weaker Cm line than the other early WC stars and therefore we were able to 
obtain a sufficiently accurate rectified line profile. Note that this line is slightly affected by 
blending with the neighbouring HellA4686 line. This is the reason why the line rises towards 
the red. However, the ClvA5016 line in that star is practically inexistent. We classify this 
St SSV as its variability level is ~ 1% of the line flux. 

For WR4, WRlll, WR5 and WR154, all WC5 and WC6 stars, no variability is de- 
tected in the weak ClvA5016 line, which is the only usable line in our wavelength range. 
We therefore classify them as stars showing no proflle variabilities (NV) but note that new 
observations of a stronger line should be obtained to verify this assertion. 

For WR135, we flnd small-scale variability only in the ClllA4650 line. The ClvA5016 
line is constant at the level of accur acy of our data set. This star has been c learly identifled 
as a typical, small-scale variable by iRobert I (Il992l ) and iLepine et al.l (Il996l ). We reach the 
same conclusion here; we classify it SSV. 

Finally, the three WC9d stars, WR106, WR119 and WR121, all reach a high level of 
variability. However, as indicated by their spectral types, these WC stars are thought to 
produce carbon-based dust in their atmosphere. So me WC dust makers have been shown 



to be binaries (e.g. the well-known system WR140; iHackwell et al.lll976l : IWilliams et al. 



19781 ). In that case, it is believed that the dust is produced within the cone-shaped shock 
formed as a consequence of the collision between the two stellar winds. As the gas flows along 
the shock cone, it eventually reaches distances where the density compression is still adequate 
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for dust formation and the intensity of the ultraviolet flux from the star is sufficiently low not 
to destroy it. For single stars, however, a definite mechanism to form dust in the wind has 
yet not been identified. Shocks with in the wind it self caused by radiative instabilities have 
been claimed to be possible sources JZubkolll993al lb[). Since all three of our WC9d stars are 
variable with a large amplitude, we believe that caution must be used, since colliding-wind 
binaries have been shown to displa y line-profile va riability in the optical with characteristics 
very similar to that of CIRs (e.g. iHill et al.l 120001 ). It is therefore possible that the WCQd 
stars in our sample are yet un-identified binaries. We plan to investigate this further by 
carrying out intensive spectroscopic and photometric monitoring campaigns of these stars. 
Only a detailed study can reveal if the changes are due to a colliding-wind shock-cone or 
to CIRs. Therefore, we classify these stars as LSV but create a separate class, the WC-late 
dust makers. 



Discussion & Conclusion 



The results of our search for optical spectroscopic variability among presumably single 
WR stars can be summarised as follows. For eight stars, WR3, WR152, WR128, WR136, 
WR131, WR158, WR108 and WR135, variability was found but at a small level compared 
to the flux of the line. This type of variability is most-likely associated with inhomogeneities 
in the dense WR wind, which move with it and reveal themselves as small-scale structures 
superposed on the broad emission lines. In more detailed studies, they hav e been found 



to mo ve fr om line center to th e edges of the lines as described for example by iMoffat et al. 



( 1l988l ) and lLepine et al.l (119961 ). As the dri ving mechanism for t hese radiatively-driven winds 
has been found to be inherently unstable (lOwocki et al.lll988l ). such variability is expected 
for all WR stars. 

Nevertheless, for 5 other stars in our sample, namely WR2, WR4, WRlll, WR5, 
WR 154 we found no significant line-profile variability. Note however that in the last 4 cases, 
our conclusion is based solely on a single relatively weak line. In fact, WRlll, is a star that 
is known to display small-scale spe ctral variability due to the presence of inhomogeneities 
in its wind and has been studied by lLepine fc MoffatI (119991 ). In any case, identifying a star 
as non-variable always depends on the quality of the data in hand. Down to what level of 
accuracy these stars are constant is, indeed, an important question. One would like to verify 
if all WR stars have inhomogeneities in their winds, if some stars present a different level of 
variability and if so, what is the cause of this difference. This kind of information can help 
us better understand the physical mechanisms responsible for these unstable outflows. 



Nine stars in our sample presented large-scale proflle variability. Of these, we have 
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identified three new candidates for CIR-type variability: WR 1, WR 115 and WR 120. Follow- 
up observations with a high temporal resolution will be required to determine if these changes 
are periodic and if so identify the period. We have done so already for WRl and present 



our results in an upcoming paper (IChene &: St-Louisll2009bl ). 



Although the results presented here are not yet for our complete sample, it is interesting 
to compare our findings to what is known for OB stars. Although the presence of DACs 



in 0-star P Cygni profiles is found to be common (e.g. iHowarth fc Prinjal 119891) . cyclical 



variability in emission lines is not found to show the same ubiquitous behavior. iMorel et al. 



( 120041 ) present a long-term monitoring campaign of the Ha line of 22 bright OB supergiants. 
Variability was found in all cases with two possible origins: either the underlying photo- 
spheric profile and/or a wind component. Cyclical changes which very likely originate from 
large-scale structures in the wind were found for only two stars in their sample (~ 14% of 
the stars). In this study we have identified 4 CIR-type variables (3 of which are new) out of 
a sample of 25, or ~ 16% , a very similar fraction. Alternatively, we can exclude the WN 8 
and WC 9 stars from this estimate; if we do so we also obtain a similar fraction, 4/19~21%. 
This striking difference between variability in P Cygni absorption components and emission 
lines is most likely due to the fact in the former we are sampling a much more physically 
restricted volume of the wind (i.e. that in front of the star) compared to the entire wind 
which can be subject to cancelling effects in parts of the wind which are physically distinct 
but that are travelling at a similar projected velocity. 

In addition to the stars showing variability thougt to be associated with CIRs, we have 
defined two distinct classes of LSV: WN8 stars which are already known to be a pecuhar 
class of WR stars and WC9d stars. In the latter case, it remains to be checked if the changes 
originate in the wind of a single star because of the presence of CIRs or if these stars are in fact 
colliding- wind binaries. If of binary origin, a very similar variability pattern is expected with 
some differences that can be identified with a detailed monitoring campaign. For example, 
strong epoch-dependancy in spite of the presence of periodic changes would suggest a single- 
star origin. Also, shock-cones translate into a very distinct and easily recognisable variability 
pattern on spectral lines, translating to at most, two large excess emissions moving around 
periodically on the emission line. 

Although our survey has revealed a relatively small fraction (~12%) of stars showing 
line-profile variability thought to be associated with CIRs among WR stars, it still provides 
some new candidates for potential estimates of WR rotation velocities. It also provides 
us with the opportunity to study in more detail the CIR phenomenon for which still very 
little is known by providing a larger pool of candidates to observe. It is clear that better 
observational constraints of this phenomenon are crucial to reach a greater understanding of 
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the physical mecanisms from which they originate. 
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